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We demonstrate the formation of bright solitons in coupled defocusing nonlinear Schrodinger 
(NLS) equation supported by attractive coupling. As an application we use a time-dependent dy- 
namical mean-field model to study the formation of stable bright solitons in two-component repulsive 
Bose-Einstein condensates (BECs) supported by interspecies attraction in a quasi one-dimensional 
geometry. When all interactions are repulsive, there cannot be bright solitons. However, bright 
solitons can be formed in two-component repulsive BECs for a sufficiently attractive interspecies 
interaction, which induces an attractive effective interaction among bosons of same type. 
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Recently, there have been successful observation [J Q 
and associated theoretical studies of two-component 
Bose-Einstein condensates (BEC). There have also been 
experimental observation ,4] of bright solitons in a BEC 
formed due to atomic attraction and related theoretical 
investigations In fiber optics true one-dimensional 
solitons are formed in the nonlinear Schrodinger (NLS) 
equation 0,0]. The solitons of BEC 4] are formed in a 
quasi-one-dimensional geometry achieved by employing 
a strong transverse trap. In either case, no solitons can 
be formed for repulsive interactions. 

In this Letter we suggest the possibility of the forma- 
tion of stable bright solitons in two-component BECs in 
the presence of repulsive interaction among like atoms 
and attractive interaction among atoms of different 
types. We find that a sufficiently strong interspecies at- 
traction can induce an effective attraction among like 
atoms responsible for the formation of bright solitons 
in two-component BECs. We also consider similar two- 
component BECs in the presence of a periodic optical- 
lattice potential, which are now routinely used in ex- 
periments on BECs 8]. We consider the formation of 
these solitons using a coupled time-dependent mean-field 
Gross-Pitaevskii (CP) equation 0. 

Because of two components such solitons of the cou- 
pled one-dimensional NLS equation are often termed 
vector solitons in nonlinear fiber optics, where two 
orthogonally-polarized pulses of different widths and 
peak powers in general propagate undistorted. The vec- 
tor solitons considered so far in nonlinear optics were 
always created in focusing (attractive) medium. To the 
best of our knowledge the present study is the first to 
consider the possibility of creation of vector solitons in 
defocusing (repulsive) medium. 

The experimental study of bright solitons in quasi-one- 
dimensional attractive systems is quite delicate due to 



the possibility of collapse in such systems, although a 
true one-dimensional system does not exhibit collapse [jj . 
The two-component repulsive BECs with interspecies at- 
traction are better suited for studying solitons as such 
systems may not easily collapse [lC} and one can have a 
controlled study of solitons. Experimentally, this could 
be realized by forming a coupled repulsive BEC in a cigar- 
shaped geometry and then transforming the interspecies 
repulsion to attraction via a Feshbach resonance 111 and 
eventually removing the axial trap so that the BEC com- 
ponents attain mobility in the axial direction like soliton 
under the action of radial trapping alone. 

Bright solitons are really eigenfunctions of the one- 
dimensional NLS equation. However, the experimental 
realization of bright solitons in trapped attractive cigar- 
shaped BECs has been possible under strong transverse 
binding which, in the case of weak or no axial binding, 
simulates the ideal one-dimensional situation for the for- 
mation of bright solitons. The dimensionless NLS equa- 
tion in the attractive or focusing case |(| 

iu t + u xx + \u\ 2 u = 0, (1) 

sustains the following bright soliton [||: 

u(x, t) — V%w sech[y/w(x — 6 + 2vt)\ 

x exp[— iv(x — 6) + i(w — vv)t + ia], (2) 

with four parameters. The parameter w represents the 
amplitude as well as pulse width, v represents velocity, 
the parameters d and a are phase constants. The bright 
soliton profile is easily recognized for v = 5 = when Eq. 
(JTJ leads to = \/2w sech(x^/w). In the case of 

the coupled focusing NLS equations: 

iu t +u xx + (\u\ 2 + \v\ 2 )u = a, (3) 
iv t + v xx + (\v\ 2 + \u\ 2 )v = 0, (4) 
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one could have the following bright vector solitons [3, 
IT^ : u{x,t) = cos(0)\/2i(;sech(a:- v /ii;) exp(it), v(x, t) = 



2 



sin(0) v 2uisech(xy / u>) exp(zt), where is an arbitrary an- 
gle. Equations and (@J are incoherently coupled as the 
coupling depends only on the intensities and is therefore 
phase insensitive. In this Letter we consider the vector 
solitons in Eqs. © and @ with defocusing diagonal 
nonlinearities. 

The time-dependent Bose-Einstein condensate wave 
function ^(r, t) at position r and time t may be described 
by the following mean-field nonlinear GP equation [|| 



lfL dl 



h 2 v 2 

2m 



+ V(r) + gn 



#(r,i) = 0, 



(5) 



with normalization J dr\^(r, t)\ 2 — 1. Here n = \~*&(r,t)\ 2 
is the boson probability density, g — 4:nh 2 aN/m, with a 
the boson-boson scattering length, m the mass and N 
the number of bosonic atoms in the condensate. The 
trap potential with axial symmetry may be written as 
V(r) — ^moj 2 (p 2 + v 2 z 2 ) where ui and vu> are the angular 
frequencies in the radial (p) and axial (z) directions with 
v the anisotropy parameter, which will be taken to be 
for axially free solitons in the following. 

In the presence of two types of bosons each of mass 
m, Eq. O ge ts changed to the following set of coupled 
equations jlfij : 



d h 2 \7 2 
in— + V(r) + gnni + gi 2 n 2 

d h 2 v 2 



* 1 (r,*)=0, 



(6) 



ih— 



dt 2m 



V(r) +52i«i +522n 2 



*a(r,t) =0. 



(7) 



Here n t = |*i(r,t)| 2 , g tj = 4nrh a^Ni/m, and i,j = 
1,2 represent the two types of bosons, and where Ni is 
the number of boson 1 and N 2 that of boson 2, ay is the 
scattering length for a boson of type i and one of type j . 

For the study of bright solitons we shall reduce Eqs. 
© and to a minimal one-dimensional form in the 
cigar-shaped geometry with v = 0. This is achieved by 
considering solutions of the type ^(r, t) — (f>i(z,t)^°'(p) 
where 



|^°>(p)| s 



7T?i 



exp 



mujp 



(8) 



The expression (JSJ corresponds to the ground state wave 
function in the absence of nonlinear interactions and sat- 
isfies 



the following one-dimensional equations |13| : 



d h d 2 

lh m + 2^d^- hM 



d h 2 d 2 



dt 
where 



2m dz 



hi\H 2 



f22\h\ 2 



0i («,*) =0,(10) 



fa{z,t) =0,(11) 



fij = 27T5i 



2?r?i' 



(12) 



In Eqs. I|10|) and (|ll|l the normalization is given by 

njM*.*)i a <fe=i- 

For calculational purpose it is convenient to reduce the 
set ()10|) and (| 1 1 1) to dimensionless form by introducing 
convenient dimensionless variables. In Eqs. (|10|l and 
(|ll(l we consider the dimensionless variables r = tuj/2, 
y = z/l, Xi = Vl4>i j with I = y/h/ [um), so that we have 
the following coupled NLS equations 



^1 

dy 2 



nn \Xi\ 



nyi |X2 1 Xi(y,T) = 0, (13) 



dy 2 



where 



AciijNi/l, i,j 



X2 (y,r)=0, (14) 



1,2. In Eqs. O 



and flU , the normalization condition is given by 
n o |Xi(l/,r)| 2 d tf = l. 

Equations (|13|l and l|14|l represent the one-dimensional 
limit of the three-dimensional equation. For solitons we 
finally have to take the nonlincarity coefficients (i j) 
to be negative corresponding to attraction. These equa- 
tions have analytic solutions only under special condi- 
tions. First, when n\ 2 = n 2 \ — 0, they become two un- 
coupled NLS equations which allow the following trivial 
soliton solutions for negative nu @: 



Xi{y,T) = v / (2wi/|nii|)sech(y x /M^exp(iT), (15) 

i = 1,2. To satisfy the normalization condition one 
should have w% = nfj/16. In Eq. 1)15(1 and below the 
parameter w is to be adjusted so as to satisfy the nor- 
malization condition. When all the nonlinear interac- 
tions are attractive (negative) and nu = n 22 = —n and 
"12 = n 2 i = —an one has the solutions [|| 



^2 1 

—V 2 ^ \p) + -muj 2 p 2 ^°\p) = h^\p),{%) Xi ( J/) T ) = y/[2w/(na + n)]sech( 2/v ^) exp(zr), (16) 



2m 

with normalization 2tt \ip(°\p)\ 2 pdp = 1. Now the dy- 
namics is carried by cj>i(z, t) and the radial dependence is 
frozen in the ground state rp(°\p). 

Averaging over the radial mode, i.e., multiplying Eqs. 
© and Q by ip^*{p) and integrating over p, we obtain 



i = 1,2, with w = (na + n) 2 /16. Also, when nu = 
0, they have the following soliton solutions for negative 
nonlinearities 



Xi{v,t) = yj '(2w i /\n Jl \)sech(yy/w~) exp(zr), (17) 
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FIG. 1: (a) The stationary functions \4>i(z,i)\,i — 1,2, for 
bosonic bright solitons vs. z for Ni = N2 = 7500, an = 1 
nm, ai2 = — 1 nm, and 022 =0.1 nm, harmonic oscillator 
length / ~ 1 /im and v = 0. The nonlinearity parameters are 
nn = 30, ni2 = W2i = —30, and 7122 = 3. 
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FIG. 2: (a) The stationary functions \4>i{z,t)\,i = 1,2, for 
bosonic bright solitons vs. z in the presence of the optical 
lattice potential V(z) = 100 sin 2 (4z) in Eqs. and (ITU . 

The nonlinearity parameters are nn = 35, 7112 ~ n%\ = —30, 
and W22 = 3. 



* 7^ h h3 = 1,2 with Wi — n^/16. This case has the 
possibility of forming the soliton due to interspecies in- 
teraction as the intra-species interaction is zero. 

It is also possible to have intra-species repulsion (posi- 
tive or defocusing no) and interspecies attraction (nega- 
tive or focusing ni2 and n2i) in order to have the fol- 
lowing soliton solutions of Eqs. (jTHjl and (jT4"|l when 
?in = 7122 = n and 7742 = 7121 = —an, (a > 1) : 

Xi(y, t) = \/[2w/(na - n)}sech(yi/w) exp(i-r), (18) 

i = 1,2 with w — (na — n) 2 /16. In this case due to 
strong interspecies attraction solitons are formed despite 
of intra-species repulsion. In all these cases the func- 
tional dependence of the two analytic solutions Xi{Ui T ) 
on y are the same. However, there are interesting numer- 
ical solutions to Eqs. (|13|) and l|14l) where the functional 
dependence of x»(2A T ), i = 1> 2, could be different, which 
we study next. Such vector solitons will have different 
widths and peak powers. 

We solve the coupled mean-field-fiydrodynamic equa- 
tions l|13[) and (|14fl for bright solitons numerically using a 
time-iteration method based on the Crank-Nicholson dis- 
cretization scheme elaborated in Ref. [l4( . We discretize 
the mean-field-hydrodynamic equation using time step 
0.0005 and space step 0.025. 

We performed the time evolution of the set of equations 
(|13fl and (|14|) introducing harmonic oscillator potentials 
y 2 in these equations and starting with the eigenfunc- 
tion of the linear harmonic oscillator problem with the 
nonlinear terms set equal to zero: Xi(y, r ) = X2(y, T ) — 
7T" 1 / 4 exp(— y 2 /2) exp(ir). During the course of time evo- 
lution the nonlinear terms are switched on very slowly 
and resultant solution iterated until convergence was ob- 
tained. Then the time evolution is continued and the har- 
monic oscillator potential terms (y 2 ) are slowly switched 



off and the resultant solution iterated 100 000 times for 
convergence. If converged solutions are obtained, they 
correspond to the required bright solutions. In the nu- 
merical investigation we take u = 2n x 100 Hz, and tub 
as the mass of 87 Rb. Consequently, the unit of length 
/ w 1 /im and unit of time 2/w ~ 3 ms. 

First we solve Eqs. (HJ) and ^ with TVi = 7V 2 = 
7500, an = 1 nm, ai2 = —1 nm and 022 = 0.1 nm. With 
these parameters the nonlinearities in Eqs. (|13fl and l|14|) 
are nn = 30, 77,12 = —30, 7121 = —30, and 77,22 = 3. 

The converged bright solitons are plotted in Fig. 1 . In 
this case the function <f>i of the first component extends 
over a longer region in space compared to the function </>2 
of the second component. It is possible to have solitons 
with different extensions in space by varying the param- 
eters of the system. The scattering length can be manip- 
ulated in the bosonic systems near a Feshbach resonance 
1 1 1 1] by varying a background magnetic field. By varying 
the scattering length and the number of atoms we could 
arrive at different values of nonlinearity parameters from 
those in Fig. 1 and thus have different extensions of the 
condensates in space. 

In the situation presented in Fig. 1, in Eq. (|13fl gov- 
erning the dynamics of boson 1 the two nonlinearities are 
nn = 30 and 7112 = —30, which may superficially indi- 
cate an effective nonlinearity of 0. However, the effective 
nonlinearity in this equation is — nn \xi | 2 — «i2|x2| 2 - Due 
to a more strongly bound soliton of boson 2 this effective 
nonlinearity could become attractive and bind the soliton 
of type 1 . 

Next we consider a two-component soliton formed in 
the optical-lattice potential V{z) = Vbsin 2 (4z) intro- 
duced in Eqs. I|13|) an( l l|14|) . In our numerical calcu- 
lation we take Vq — 100, nn = 35, ni2 = n2i = —30, 
and n22 = 3. In this case in the wave function of boson 
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FIG. 3: The function [0j(2,t)| for bosonic bright solitons vs. 
z and t for (a) boson 1 and (b) boson 2. At t = Ni = 
N2 = 7500, an = 1 nm, 012 = —1.7 nm, and a22 = 0.1 
nm, harmonic oscillator length I ~ 1 (jm and v = 0. The 
nonlinearity parameters at i = are tin = 30, nvi = TI21 = 
—51, and 7122 =3. At t = 100 ms (marked by arrows) the two 
bright solitons are set into breathing oscillation by suddenly 
jumping the nonlinearities ni2 and 7121 to n.12 = TI21 = —45. 



1 prominent wiggles are formed due to the optical-lattice 
potential. As the spacing of the optical-lattice sites are 
relatively large no wiggles are formed in the more local- 
ized wave function of boson 2. Equation (|13|) taken sep- 
arately has nonlinearities «n = 35 and nn = —30 corre- 
sponding to an apparent overall repulsion. However, the 
strongly bound soliton of boson 2 and the interspecies 
attraction as well as the periodic optical-lattice potential 
aid in the formation of the soliton of boson 1 . In this con- 
nection it should be noted that the optical-lattice poten- 
tial aids in binding the soliton in the presence of nonlin- 
ear interspecies attraction. The optical-lattice potential 
alone cannot bind a soliton in the absence of nonlinear 
attraction both in one- and two-component BECs. 

Finally, we consider the stability of these solitons under 
a small perturbation. For this purpose we consider the 
solitons of Eqs. (|13fl and l|14|) formed for nonlinearities 
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FIG. 4: Root mean square size (z)rms vs. time of the two 
solitons of Fig. 3 set into breathing oscillation by suddenly 
jumping the nonlinearities at t = 100 ms. 



nn = —30, ni2 = n<2.\ = —51, and 7122 = 3. The wave 
functions 4>i{z, t) of these solitons are shown in Figs. 3 (a) 
and (b) for i = 1,2, respectively. To test their stability 
under small perturbation, after their formation, the non- 
linearities tt-12 and 71,21 are suddenly changed to —45 at 
t = 100 ms so that the solitons are set into motion. Such 
change in the nonlinearity can be achieved by a jump 
in the scattering length by manipulating a background 
magnetic field near a Feshbach resonance 11]. The soli- 
tons are found to execute stable non-periodic breathing 
oscillations. The stability of these solitons after the per- 
turbation is applied is demonstrated in Figs. 3. 

To further illustrate the stability of the oscillation of 
the solitons of Figs. 3 we plot in Fig. 4 the root mean 
square (rms) size (z) rms vs. time. Sustained oscillation 
for a very long time illustrates the stability of the soli- 
tons. From Fig. 4 we find that for t < 100 ms the rms 
sizes of the two solitons are constant. However, after the 
application of the repulsive impulsive force at t = 100 
ms by reducing the attractive nonlinearity, the rms sizes 
suddenly jump to a larger value and execute stable os- 
cillatory dynamics. This stable oscillation guaratees the 
stationary nature of the solitons under small perturba- 
tion. 

The present investigation has consequences in generat- 
ing bright vector solitons in directional couplers of nonlin- 
ear fiber optics ,7]. Vector solitons are solutions of cou- 
pled one-dimensional NLS equations with the property 
that the orthogonally polarized components propagate 
in a birefringent fiber without change in shape. In vector 
solitons an input pulse maintains not only its intensity 
profile but also its state of polarization even when it is 
not launched along one of the principal axes of the fiber. 
We have investigated the possibility that two orthogo- 
nally polarized pulses of different widths and different 
peak powers propagate undistorted in birefringent fibers. 
The present investigation suggests that it is possible to 
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have bright vector solitons in the coupled NLS equation 
with diagonal defocusing (repulsive) nonlinearity and off- 
diagonal focusing (attractive) nonlinearity, which should 
be of interest in nonlinear fiber optics. 

We use a coupled set of time-dependent mean-field GP 
equations for a two-component repulsive BEC to demon- 
strate the formation of bright solitons due to interspecies 
attraction. The interspecies attraction can neutralize the 
intra-species repulsion and induce an effective attraction 
in the mean-field GP equations responsible for the forma- 
tion of bright solitons. An attractive interspecies interac- 
tion is necessary for the formation of the bright solitons 
as the diagonal nonlinearity na in the mean-field equa- 
tions is taken to be repulsive (positive). In mean-field 
equations 1)13(1 and (f 1 4|) nus are positive (repulsive) and 
riijS are negative (attractive), so that the overall contri- 
bution of the nonlinear terms in these equations become 
attractive to support the bright solitons. We have also 
established the formation of these solitons in the presence 
of a periodic optical-lattice potential in an entirely dif- 
ferent shape and trapping condition from a conventional 
soliton. 

In view of the present study the appearance of bright 



solitons in multi-component repulsive BECs seems pos- 
sible in quasi-one-dimcnsional geometry. Bright solitons 
have been created experimentally in attractive BECs in 
three dimensions in the presence of radial trapping only 
without any axial trapping Also, there have been 
experimental studies of multicomponent BECs 0, Q]. 
Hence, bright solitons can be created and studied in the 
laboratory in the presence of radial trapping only in a 
two-component repulsive BEC supported by interspecies 
attraction and the prediction of the present study veri- 
fied. We have also suggested the possibility of the for- 
mation of similar coupled solitons in fiber optics in one 
dimension. 

Note added in proof: After the completion of this in- 
vestigation we have known about another similar study 
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